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ABSTRACT 

Two parametric, ensemble cloud models have been tested 
with data obtained during GATE (1974). The first model 
is an adaption of the Arakawa-Schubert scheme which uses 
an entraining jet to represent individual cumulus cloud 
types. The second model consists of an ensemble of cylin- 
drical cells to represent the convective cloud field. 

IMTODUCTION 

An important feature of synoptic- scale numerical prediction 
models and for planetary- scale general circulation models is the 
parameterization of cumulus convection. Fbr integrations over 
short time periods, it is possible to treat subgrid-scale con- 
vection in an extremely artificial way, or ignore it altogether. 
However, when the period is large, and especially in the tropics, 
convection is an important mechanism for the vertical transport 
of heat, moisture and momentum. Even when the period is short 
(<1 day) and the large-scale flow is not influenced by con- 
vective activity, the prediction of convective rainfall and 
cloudiness are still important objectives of the local weather 
forecast and climate of a region. 

In preparation for a study of parameterized convection in 
a predictive model where physical processes are not free to be 
chosen arbitrarily, we have studied several parameterization 
models in a diagnostic sense, using observations of synoptic-scale 
variables to determine the characteristics of convection rather 
than the reverse. If these models give reasonable results, 
their ultimate application will be in a prediction scheme in 
NWP and GCM models. 

The first part of the study investigates the sensitivity of 
the Arakawa-Schubert parameterization method. Lord (1978) has 
completed a more extensive study using a different numerical 
procedure ("overadjustment"). We have also compared our more 
limited results with those obtained using an overadjustment 
procedure by Schubert and Silva-Diaz (personal communication) . 
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The second! part of the study is a test of an ensemble of 
cellular convective elements (Cheng, 1978; Rodenhuis and Cheng, 
1979). The results of this model are compared to those of 
Nitta (1977) who also used the entraining jet cloud model of 
Arakawa and Schubert (1974) to explain the residual term in a 
budget calculation. 

RESULTS 

Sensitivity Tests of the Entrainment Model 

Figure 1 shows the results of computing the detrainment 
level for each cloud type vrtiich is identified by a constant 
entrainment coefficient, The input data is the composite 
of Nitta (1977) for unorganized convection in GATE (6 September 
1974). (The radiational heating rate, Q R = 0 for these studies.) 
Three cases are shown corresponding to different conditions at 
the cloud base: 



Fig. 1. The detrainment level of model clouds as a function of 
the detrainment rate (X) . Case 1 (solid) ; Case 2 
(dashed) ; Case 3 (dash-dotted) . 
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Case Is Cloud base and top of the mixed layer all 
identical; cloud elements are neutrally 
buoyant at cloud base. 

Case 2: Cloud base is identical with Case 1, but the 

top of the mixed layer is assumed to be 25 mb 
below? cloud elements have the properties of 
the environment within the mixed layer. 

Case 3: Cloud base and the top of the mixed layer are 

identical; cloud elements have the same properties 
as the environment at the cloud base (Ah=As y =0) . 

These rather small changes in conditions which are applied 
at the cloud base strongly influence the spectra of convective 
cloudiness as shown in Fig. 2. Only in Case 1 is there a large 
number of small cumuli which are considered to be so essential 
for the moisture budget of an atmosphere with imbedded deep 
convection. 



Fig. 2. The cloud mass flux spectrum as a function of detrain- 
ment level for the three cases shown in Fig. 1. 


Figure 3 shows a comparison of our results (University of 
Maryland [UM] ) with those from a calculation of W. Schubert and 
P. Silva-Diaz (Colorado State University [CSU] ? personal communi- 
cation). For this comparison, identical input data was used. 

The differences are due almost entirely to the method of inte- 
gration in the calculation of A(z): 

h* = h (z,X) = [h (z^,X) +xf n(z",X)h(z')dz'] (1) 

C n(Z,A) C B / 
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where h* and h are the moist static energy in the environment 
C and cloud, respectively. The * indicates 
scattered values. 

n(z,A) is specified exponential function of height 

and A. 

z is the height of the cloud base. 

B 

Note that there is a gap in the spectrum of the CSU 
results between 875 and 925 mb (dashed line) which does not 
occur in the UM calculation. 



Fig. 3. The detrainment level of model clouds as a function of 

detrainment rate (A) computed with two different methods. 


This discrepancy accounts for some of the differences in 
the cloud spectra seen in Fig. 4. In addition, the UM results 
have been simply minimizing the residual term in the equation 
(the equilibrium hypothesis): 


where 


and 


dA 

dt 



dA + F = 0 


A is the cloud work function 
K is the cloud :cloud interaction matrix 
F represents the large-scale forcing function 
M b is the cloud base mass flux. 


( 2 ) 


However, Schubert uses_a different method (described by lord, 
1978), which modifies F to accommodate changes in A which 
occur over a small arbitrary period of time. At. Thus, 
Schubert's results depend on At, as can be seen in Fig. 4 
for the 2 cases of CSU results. 


166 




Fig. 4. The cloud mass flux spectrum as a function of detrain- 
ment level for three cases. 

Note that the UM results do not prohibit low cloud types 
(Fig. 3) , but they are not required for the solution of (2). 
However, we do require high clouds detraining at 175 mb. By 
way of control, the CSU results of 875 mb (dotted curve) prohibit 
cloud tops at that level, but not more shallow clouds. Therefore, 
a large number of shallow clouds are required to satisfy the 
equilibrium condition, at least for At=10 min. 

Comparison with the Cellular Model 

Figure 5 (top) shows the results of cloud spectra obtained 
by using an ensemble of cellular convective elements (Rodenhuis 
and Cheng, 1979). This result also differs from those in Figs. 2 
and 4, because an entirely different method is used to compute 
the cloud number density and mass flux. Rather than imposing 
the equilibrium hypothesis (2), the net-heating-and-^noisture- 
deficit budget is satisfied: 


°1-W -3p (e V- l |j ff,) < 3 > 

where = net cumulus heating 

Q 2 = net cumulus moisture deficit 
Q r = radiational heating rate 
F 0 = net cumulus heat flux 
Fg = net cumulus moisture flux. 
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Fig. 5. A comparison of cloud base mass flux determined frcm 

the cellular model (upper) and determined by Nitta (1977) . 
Input data (6 September 1974, GATE) for both results are 
identical . 

Although there are some similarities between Figs. 2, 4, 
and 5, there are serious discrepancies in every case. (Sane 
differences are due to different values of Q^. ) 

Comparison with the Entrainment Model 



CLOUD HEIGHT, km 


The entrainment cloud model has also been used with (3) 
by Nitta (1977) to determine the cloud spectra. These results 
are shown in Fig. 5 (lower). The relatively large number of 
lower cloud agrees (perhaps fortuitously) with Case 1 of Fig. 2. 
There is some similarity between the upper and lower parts of 
Fig. 5, but the bimodal distribution does not occur with the 
cellular model because of the increased mass flux of middle-level 
clouds. 

CONCLUSIONS 

The spectra of the cloud base mass flux has been calculated 
by two different methods: (1) by an application of the equilib- 

rium hypothesis, equation (2), and (2) by using the bulk budget 
method, equation (3). In each cse we have shown two independent 
calculations using different computational schemes (UM and CSU), 
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or using different models (entrainment model and the cellular 
model). Our results using the equilibrium hypothesis show that 
the spectra are very sensitive to the conditions at the cloud 
base. Furthermore, the results depend on the computational 
scheme which is used. Finally, the spectra do not agree very 
well with those calculated from the budget method, although 
differences in Q R account for some of this discrepancy. 

When the cellular model is used with the budget method, 
the results are in general agreement with Nitta's spectra, 
although there is not a strong bimodal mass flux distribution. 

The rainfall estimate with the model is 10.21 mm/day, as compared 
to 10 mm/day from the average of the shipboard observations 
during the period. Other evidence (Cheng, 1978) shows that the 
cellular model behaves reasonably over a wide range of conditions. 

In addition, the ensemble cellular cloud model can diagnose 
the cloud areal fraction and number density of each cloud type. 
Figure 6 shows these new spectra for the same period of analysis 
as has been used throughout this study. These results may be 
used to compare with independent, direct measurements of cloud 
dimensions from radar or satellite radiometers. This is an 
extremely useful link between cloud parameterization in numerical 
models and direct observations of clouds, especially from satellite 
platforms. 
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Fig. 6. Spectra of cloud areal fraction and number density as a 
function of closed height determined from the ensemble 
cellular model. 
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